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ABSTRACT 

The submillimcter opacity of dust in the diffuse interstellar medium in the Galactic plane has been 
quantified using a pixel-by-pixel correlation of images of continuum emission with a proxy for column 
density. We used multi-wavelength continuum data: three BLAST bands at 250, 350, and 500 [im 
and one IRAS at 100 /im. The proxy is the near-infrared color excess, E(J — K s ), obtained from 
2MASS. Based on observations of stars, we show how well this color excess is correlated with the 
total hydrogen column density for regions of moderate extinction. The ratio of emission to column 
density, the emissivity, is then known from the correlations, as a function of frequency. The spectral 
distribution of this emissivity can be fit by a modified blackbody, whence the characteristic dust 
temperature T and the desired opacity er e (1200) at 1200 GHz or 250 /Ltm can be obtained. We have 
analyzed 14 regions near the Galactic plane toward the Vela molecular cloud, mostly selected to 
avoid regions of high column density (iVn > 10 22 cm~ 2 ) and small enough to ensure a uniform dust 
temperature. We find <7 e (1200) is typically 2 to 4 x 10~ 25 cm 2 H _1 and thus about 2 to 4 times larger 
than the average value in the local high Galactic latitude diffuse atomic interstellar medium. This 
is strong evidence for grain evolution. There is a range in total power per H nucleon absorbed (and 
re-radiated) by the dust, reflecting changes in the strength of the interstellar radiation field and/or 
the dust absorption opacity. These changes in emission opacity and power affect the equilibrium 
T, which is typically 15 K, colder than at high latitudes. Our analysis extends, to higher opacity 
and lower temperature, the trend of increasing a c (1200) with decreasing T that was found at high 
latitudes. The recognition of changes in the emission opacity raises a cautionary flag because all 
column densities deduced from dust emission maps, and the masses of compact structures within 
them, depend inversely on the value adopted. 

Subject headings: Balloons - dust, extinction - evolution - Infrared: ISM - ISM: structure - Submil- 
limeter: ISM 
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1. INTRODUCTION 

Observations of submillimeter dust emission are a 
prime means for determining masses in the interstel- 
lar medium, including in compact clumps and cores of 
star-forming regions. However, there remains consider- 
able systematic uncertainty because the opacity er e , the 
dust emission cross-section per H nucleon, is not well 
known. The value of a c is best determined for the dif- 
fuse atomic interstellar medium at high Galactic lati- 
tude ( Boulang er" et al.lll99"5 iPlanck Collaboration et all 
l2011al)~ The main goal of this paper is to quantify a e 
empirically in different environments near the Galactic 
plane, characterized by higher column and spatial den- 
sity and at least in part molecular. 

Detailed models of interstellar grains dDwek et al.l 
fT997tlLi fc Dram^l2lfflirCommegne etalMmF are best 
constrained in the local diffuse and largely atomic inter- 
stellar medium. The value of o~ e can be calculated for a 
given model of dust grains, depending upon the compo- 
sition and grain structure but fortunately not so strongly 
on size or shape. This calculation involves the emission 
cross section per gram of dust, k v , and the dust-to-gas 
mass ratio, r. Thus these models are further constrained 
by and/or checked for the consistency of the product rn v 
with the empirical value of a c (Equation (|4])) at high 
Galactic latitude. 

The value of the submillimeter opacity is likely to 
change with environment, through differences in com- 
position, structure, and even dust temperature. Cer- 
tainly there are changes in the optical-ultraviolet extinc- 
tion curve that indi cate and even quanti fy some aspects 
of dust evolution (Kim & Mar tini Tl996l and references 
therein). Pioneering work on NGC 7023, correlating sub- 
millimeter emission with a measure of ultraviolet extinc- 
tion, demonstrated empirically that the opacity might 
be higher by a fac tor of two in that denser environment 
( Hildc hrandl 1 1983T) . However, the uncertainty was much 
too large (a factor three or four) to be definitive. 

Where, when, and how evolution happens are un- 
known, but it seems plausible that in denser and colder 
regions, where dust grains acquire ice mantles, the dust 
is more susceptible to aggregation on a relevant time 
scale. Theoretical dust models can be used to explore 
the complexities of the observable consequences (extinc- 
tion, emission) of dust evolution i n denser environm ents 
(e.g.. lOssenkopf fc Henn"n7elll994l : lOrmel et al.ll2011b see 
also Section 18. 3p but the evolution is undoubtedly too 
complex for precise environment-specific predictions ab 
initio. Thus we take an empirical approach to quan- 
tifying the opacity in different environments of higher 
column and spatial density. 

The basis of our study is a careful correlation of the 
emission in submillimeter dust continuum images ob- 
tained by the Balloon-borne Large Aperture Submil- 
limeter Telescope (BLAST), together with images from 
IRAS, with a proxy for column density, the near- infrared 
color excess that is obtained from analysis of 2MASS3 
data. The slopes of these correlations are emissivities. 

18 The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 



BLAST mapped sim ultaneously at 250, 350, and 500 fim 
(jPascale et al.ll2008D and so together with IRAS 100 ^m 
data the spectral coverage is quite good in the region 
of the peak of the spectral energy distribution (SED) 
of the emissivities, allowing the characteristic dust tem- 
perature T to be constrained via a modified blackbody 
fit. Knowing T is essential for recovering the dust opac- 
ity from the dust emissivity. Accompanying goals are to 
assess the uncertainties, to acknowledge explicitly some 
"known unknowns," to show how our results extend the 
range of environments in which there is a reasonable cal- 
ibration of the opacity, and to explore trends between 
variations in T and o~ e . 

Our paper is organized deliberately to isolate the steps 
in our approach and to identify potential sources of sys- 
tematic error. In Section [2] we briefly discuss mecha- 
nisms of diffuse emission and the quantitative relation- 
ships of dust emissivity and opacity to mass estimates 
from submillimeter continuum emission. We introduce 
the BLAST imaging in Section [3] and make our first es- 
timate of T from the relative SED. In Section [4] we dis- 
cuss color excess measured with 2MASS, a surrogate for 
column density. The fundamental correlation of emis- 
sion with color excess is presented in Section [SJ where 
we then use the SED of the emissivities to obtain the 
amplitude and a second measure of T, corroborating the 
first. In Section [6] we determine the ratio of hydrogen 
column density to infrared color excess, E(J — K s ), from 
stellar data. This allows us in Section [7] to determine the 
desired opacity a c from the above amplitudes and tem- 
peratures and to compute other important parameters 
like P, the power emitted (absorbed) per H nucleon. We 
assess the errors and in Appendix [5] the impact of possi- 
ble systematic changes in the spectral dependence of the 
opacity with T. We show that the parameters T, cr e , and 
P change significantly from region to region. In Section[S] 
we discuss the systematic interrelationships among these 
changes, bringing in for comparison other estimates of 
opacity drawing on a brief summary of the literature in 
Appendix [B] We comment on underlying uncertainties, 
range of applicability, and efforts at theoretical modeling. 
Finally, Section |H] gives our conclusions and anticipates 
future work. 

2. SUBMILLIMETER DUST EMISSION 

BLAST maps of thermal dust emission measure surface 
brightness I v (MJy sr _1 ) and hence, for optically-thin 
emission, the mass column density of the dust, which 
is a fraction r (the dust-to-gas mass ratio) of the total 
mass column density. Thus, more technically, 

I v = M d K u B v (T) = r^m n N n K v B v {T) = r v B u (T), (1) 

where k„ is the dust mass absorption (or emission) co- 
efficient (cm 2 g^ 1 ), often called the opacity, B V (T) is the 
Planck function for dust temperature T, /i is the mean 
molecular weight, Nn is the total H column density (H 
in any form) , and t„ is the optical depth of the column 
of material. If T changes along the line of sight, then 
Equation ([1]) has to be generalized appropriately, but 
obviously the analysis is more straightforward if strong 
T gradients are avoided, as we will do. 

For a (small) region of uniform T the variation in 
brightness in a map tracks changes in optical depth along 
the different lines of sight. In fact, as Equation ([I]) shows, 
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a map of optical depth could be obtained directly by di- 
viding the image by the Planck function B u (T) if the ap- 
propriate dust temperature can be obtained. However, 
because the "zero point" of the BLAST emission maps 
is not known we do not use this direct approach here. 
The emissivity of a column of interstellar material, 

e e (u) = I V /N H , (2) 

is an observable prerequisite to quantifying the desired 
opacity. By analogy with Equation ([2]), for a column of 
interstellar material characterized by its extinction, we 
can define 

e c (u)=I u /E(J-K s ), (3) 

where E{J — K s ) is the color excess from differential ex- 
tinction between the J and K s passbands (of 2MASS); 
e c {v) is the observable addressed below. 

The opacity <J e [y) of the interstellar material medium 
is ultimately a property of the grain material as can be 
seen from the relationships 

a c (v) = Tu/N-r = e c {v)/B v {T) = ^m^rK^. (4) 

Quantifying the opacity also requires knowledge of T, 
which can be obtained from the multi-wavelength SED 
of the emissivities as described below. 

In what follows we parameterize the spectral depen- 
dence as k„ = Ko(^/^o)^, with fiducial frequency vo = 
1200 GHz (A = 250 /iin) and f3 = 1.8, to compare 
directly with the value of <7 e (1200) for the diffus e in- 
terstellar medium given by Boulang er et al.l (|1996f ) and 
iPlanck Co llaborati on et al.l ()2011aD . 

For identifiable objects in a map (peaks of emission) 
the integral of l v over source size gives a flux density. 
When flux density is integrated over v across the entire 
SED, and the distance is known, a luminosity is obtained. 
One does not have to know the opacity to interpret this 
observable; however, both opacity and temperature are 
required to obtain the mass of the (compact) object. By 
analogy, a related quantity for diffuse emission is the 
power per H emitted by dust grains (equal to the ab- 
sorbed power), computed by integrating the emissivity 
over v. 

P = J 4ne c (is)dv = J \Tto c {v)B v (T)dv. (5) 

Again this can be computed from the observable set of 
e c (i>) values, without knowledge of the opacity. It can 
also be seen from our parameterization of the frequency 
dependence of the opacity that P varies as (T e (1200)T 4+ *. 

For t he diffuse atomic high latitude interstellar 
medium ()Planck Co llaborat ion et al.ll2011aj) , typical val- 
ues are /3 = 1.8, T = 17.9 K, <7 e (1200) ee er = 
1.0 x 10~ 25 cm 2 H" 1 (rK = 0.043 cm 2 gm" 1 ) and 
Pq = 3.8 x 10~ 31 W H -1 . Perhaps more memorable, 
expressed in solar units this power is close to unity: 
-Po/wh ~ 1.2 L q /Mq. Note that P is not solely a prop- 
erty of the dust; because the thermally-emitting dust is in 
radiative equilibrium, P depends directly on the strength 
of the interstellar radiation field in which the dust is im- 
mersed. In the high latitude analysis, the opacity is cal- 
ibrated through empirical correlation of dust emission 
with atomic hydrogen column density. It is found, sur- 
prisingly, that the opacity is not constant even in that 
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Figure 1. BLAST 250 fim dust emission map spanning the Vela 
Molecular Ridge. This is projected on the grid of the 2MASS-based 
color-excess map in Figure [2J at a degraded resolution and 1.'5 
per pixel. Outlined are several small rectangular regions, aligned 
roughly along the scans, which are analyzed independently. These 
were selected to avoid gradients in dust temperature, compact 
sources, and (molecular) regions of the highest column density. 

high latitude environment, that one would expect to be 
simple. Our new study extends such studies to regions 
that have higher column density and are at least partially 
molecular. 

3. SUBMILLIMETER OBSERVATIONS: BLAST IMAGING 
OF VELA 

Wi th BLAST (|Pascale et al.l [20081: [Truch et al.l [20081 
2009) we surveyed 50 deg 2 in Vela for 10.6 hours dur - 
ing the December 2006 flight (jNetterfield et al.l 120091) . 
BLAST06 produced diffraction-limited images with res- 
olutions (full width at half maximum: FWHM) of 36", 
42", and 60" at 250, 350, and 500 ^m, respectively 
(jTruch et a l. 2009) . For the pixel-pixel correlations with 
the color excess maps (Section [5]), we have convolved the 
BLAST images to 2.'1 resolution and regridded to 1.'5 
pixels. The 250 /im map used is shown in Figure [TJ 

The observations were performed by scanning the tele- 
scope in azimuth at a speed of 0?2 s . The combination 
of high scan speed and low 1// knee, together with the 
multiple cross-linking an d common- mode remov al in the 
map- maker SANEPIC (jPatanchon et all 120081) retains 
diffuse low spatial frequency emission (the DC level is 
removed, however). Ideally the cross-linking scans would 
have been orthogonal, but solar position constraints re- 
sulted in a small angle range, with scans being oriented 
roughly along constant declination. Consequently small 
drifts in the baseline produced a low spatial frequency 
undulation in the cross-scan direction, readily apparent 
by comparing results from different independent map- 
makers. We have examined rectangular subregions of 
sufficiently small dimension in declination to mitigate 
any potential effects of this artefact (much taller regions 
were, in fact, examined too and the results were found 
to be robust). 

The subregions examined are shown in Figures [1] and [2] 
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Figure 2. Like Figure [l] but for the column density from extinc- 
tion, the color excess E(J — K B ). The intensity scale has been 
labeled in Ay, which is perhaps more intuitive, but as discussed 
in Section [4] this conversion is only approximate and should not be 
the basis for quantitative calculations. 



Table 1 

Rectangular regions analyzed 



ID 


a 


(J2000) 


8 (J2000) 


I 


b 


Size 






h:m:s 


d:m:s 


O 


O 


/ 


X 


/ 


1 


8 


45:39.5 


-44:24:22.9 


264.061 


-0.868 


64.5 


X 


30.0 


2 


8 


48:38.6 


-41:58:41.9 


262.511 


1.078 


61.5 


X 


28.5 


3 


8 


51:00.2 


-46:02:22.4 


265.923 


-1.166 


64.5 


X 


22.5 


4 


8 


51:58.6 


-41:08:50.7 


262.264 


2.087 


76.5 


X 


27.0 


5 


8 


52:04.5 


-45:34:11.9 


265.680 


-0.722 


51.0 


X 


22.5 


6 


8 


54:42.2 


-44:55:47.2 


265.486 


0.043 


87.0 


X 


16.5 


7 


8 


59:38.4 


-46:03:41.7 


266.908 


-0.036 


51.0 


X 


24.0 


8 


9 


04:59.1 


-46:15:16.4 


267.666 


0.529 


79.5 


X 


19.5 


9 


9 


05:48.6 


-44:25:41.8 


266.407 


1.860 


61.5 


X 


30.0 


10 


9 


07:14.3 


-45:47:51.4 


267.590 


1.126 


72.0 


X 


22.5 


11 


9 


08:16.0 


-44:58:07.5 


267.101 


1.818 


94.5 


X 


24.0 


12 


9 


12:53.4 


-46:42:54.4 


268.928 


1.216 


85.5 


X 


33.0 


A 


8:55:28.1 


-46:00:21.9 


266.394 


-0.550 


25.5 


X 


16.5 


13 


9:01:12.0 


-44:53:13.3 


266.204 


0.943 


34.5 


X 


13.5 



and specified in Table Q] The typical length of a rectan- 
gle, 1°, corresponds to 12 pc at the distance of th e Vela 
Molecular Ridge, 700 pc (|Murphv fc Mavlll99lD . The 
Galactic latitude range covered by the centers of the dif- 
ferent rectangles is —1° to +2°. The maximum sepa- 
ration in longitude is about 7° (85 pc). The different 
rectangles sample distinct environments. 

To avoid areas with strong gradients in dust temper- 
ature, we concentrated on relatively diffuse regions of 
moderate brightness, without compact sources, not the 
high column (and spatial) density molecular clouds where 
conditions might be less uniform along the line of sight 
and where the uncertainties in determining the color ex- 
cess E(J — K s ) are larger. 

However, suspending our caution, we examined in 
addition two promi nent higher column den sity regions 
traced even in 13 CO (|Yamaguchi et al.lll999f) : see rectan- 



gles A and BFI lNetterfield et all (f2009h derived cr o (1200) 
for the clump in rectangle A by comparing the integrated 
BLAST emission with the mass derived from CO obser- 
vations (|Yamaguchi et al.1 119991 ) . providing an indepen- 
dent check on both approaches. 

Note that for absolute measures of column density, 
pixel by pixel, we would have to restore the zero point 
(DC le vel) of the BLAST ma ps, as we did for the Cas A 
region (jSibthorpe et al.ll2009fl . However, that is not nec- 
essary here, since we are exploiting the spatial correla- 
tions of dust emission with color excess over the individ- 
ual rectangles. A corollary is that we obtain the proper- 
ties of the dust that is producing the spatial variations 
within each rectangle; a uniform distribution would not 
be detectable. 

3.1. Dust Temperature from the Diffuse Emission 

Small scale structures are remarkably well correlated 
across the three BLAST bands. For a sufficiently large 
and homogeneous region an estimate of the characteris- 
tic temperature can be obtained via pixel-by-pixel cor- 
relations of images with respect to some reference im- 
age (here BLAST 250 H- We used the IDL routine 
SIXLIN to perform the regressions, a nd adopted the b i- 
sector slope as our specific estimator (jlsobe et al.i ri990) . 
The slopes of these correlations describe the relative SED 
of the emission that is changing in common in these im- 
ages. 

The SED for cold dust emission at temperature ~ 
15 K and emissivity index f) = 1.8 peaks at 200 /Ltm, 
and so in principle the dust temperature could be de- 
termined from the curvature in the SED through the 
three submillimeter passbands of BLAST. Nevertheless, 
in practice it is always preferable to have broader wave- 
length coverage, particularly on the short wavelength 
side of the peak. At the shorter wavelengths we have 
examined the correlation slope of the 100 /im image 
from IRA S with BLAST. We used the reproc essed IRIS 
product (jMiville-Deschenes fe Lagachd [2005). Because 
there was a hot point source in rectangle 1 and another 
in rectangle 4 we used a version of the IRIS images 
from which the sources have been removed (these hot 
sources are not prominent at the BLAST passbands). 
The IRIS 100 /im image has 1 '5 pixels in common with 
the other maps used but a lo wer resolution of about 4' 
( Mivillc-Dcsch jnes et al.ll2002ft . We considered using an- 
other IRAS prod uct, the higher resolution HIRES image 
(|Cao et al.M1997h . HIRES was orig inally designed for re- 
solving crowded compact sources (jAumann et al.lfl990T ) 
and although empirically it improves the res olution of 
diffuse emission too (e.g., iMartin et al.l H994). this has 
not been quantified at the few percent level. We have, 
however, convolved the HIRES map (resolution 2') to 
the IRIS resolution and correlated it with the original 
within the rectangles examined. The slope was typically 
0.96. We concluded that using IRIS might result in a 
systematic error of a few percent, but that this was ac- 
ceptable given the much larger calib ration error of 13% 
(jMiville-Deschenes fc Lagachd I2005D . The effect of such 

19 Rectangle B is in the main Vela C cloud and unavoidably 
contains some cold compact sources that, although peaks in column 
density, might potentially produce extra scatter in the correlation 
with the color excess map, as discussed below in Section [5] 
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systematic errors is explored in Section 17.41 

An example of the relative SED is shown in Figure 02 
for rectangle 7. Note that the fitting function implic- 
itly passes through unity at 250 /im (indicated by a dia- 
mond); there is no datum there used explicitly in the fit 
of the relative SED. 

The diffuse emission in the submillimeter and mid- 
infrared wavelengths comes from two different dust com- 
ponents, distinguished principally by their size dis- 
tribution, namel y Big Grains (BGs) and Very Small 
Grains (VSGs) (IDesert et all [TOM ILi fc Drainel EfflOTl : 
iCompiegne et al.l 12011 ). The BGs, in thermal equilib- 
rium with the ambient radiation field, account for most 
of the dust mass and therefore most of the longer wave- 
length emission. VSGs are small enough to experience 
non-equilibrium heating and so broaden the spectrum 
toward shorter wavelengths, beyond the spectral peak 
of the BG emission. VSGs comprise a relatively low 
fraction of the total dust mass, even less in dense re- 
gions, and typically their excess emission over what is 
expected from equilibrium BGs alone appears at 60 /jm 
and shorter wavelengths; in the fields studied here this 
emission is very faint. Therefore, in this exploratory work 
we adopted a single temperature SED and fit only data 
for wavelengt hs 100 11m and l onger, using the IDL rou- 
tine MPFIT (|Markwardtl2009lL Even with fixed P = 1.8, 
this grey-body functional form provides an acceptable fit 
to the data (e.g., Figure E] the best-fit equilibrium tem- 
perature for this rectangle is 15.3 ± 0.4 K). 

The observed low slope of the correlation of the 100 /im 
IRIS image with the 250 /im BLAST image constrains 
the dust temperature to be low. The degree of corre- 
lation between the 100 /im and the 250 /im emission 
is somewhat less than it is between different BLAST 
bands. This decorrelation is probably due to a range of 
grain temperatures within the volume sampled, possibly 
including a contribution from non-equilibrium emission; 
temperature changes have non-linear effects in the Wien 
tail. The correspondingly larger uncertainty in the slope 
of the correlation, by about a factor two, means that the 
100 fim datum has less weight in fitting the SED. 

4. OBSERVATIONS OF COLOR EXCESS 

To determine the dust opacity from the submillimeter 
emission we need an independent tracer of column den- 
sity, here the near-infrared color excess. This can be es- 
timated using J, H, and K s data from the 2MASS point 
source catalog, through a variety of techniques. The 
color excess map used here was derived by the "AvMAP" 
procedure which calculates the average near-infrared 
reddening of stars with a method adapted from ear- 
lier analyses (|Lada et al.l 119941 : lLombardi fe Alvesl 120011 : 
ICambresy et al.ll2002h. w i th im provements as described 
briefly in Schneider et al.1 (j2006[ ) and more completely in 
iSchneider et al.l (|2011h . For diagnostic purposes, we have 
used this procedure to produce maps of both E(J — H) 
and E(H — K s ). These maps have resolution of about 
2.'1 and 1.'5 pixels. 

We have checked these maps against those more 
recently published for the whole sky. They cor- 
relate well with the ex tinction maps produced by 
iRowles fc Froebrichl ([2009) using a median near-infrared 
color excess technique and w ith the color excess maps 
produced by iDobashil (|2011f ) using the "X" percentile 
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Figure 3. Relative SED from correlations in rectangle 7 of 
BLAST and IRIS data with BLAST 250 /im. Implicitly the value 
is unity at 250 (im (indicated by a diamond); there is no datum 
there used explicitly in the fit of the relative SED. With = 1.8 
the best fit temperature is 15.3 ± 0.4 K. 

m ethod. Fo r exam ple, our E(H — K s ) compared to that 
of IDobashil (pOlTI ) is 1.08 ± 0.01 for the moderate col- 
umn densities of interest in our study (corresponding to 
Ay < 10), but with considerable deviation beyond that, 
where extinction is harder to determine 

Our two color excess maps E(J — H) and E(H — K s ) 
are very tightly correlated, with a slope of 1.77±0.01 over 
a range up to E( J — K s ) of 5. The rectangles examined 
here had maximum E{ J — K s ) typically less than a fifth 
of this value (Tabled . Even over the whole field, there is 
no curvature in the correlation as might be diagnostic of 
a medium where dense clumps below the resolution limit 
of "AvMAP" systematically dimmed stars to beyond the 
completeness limit at J first, biasing E(J — H) to be low. 

Theoretically, the E(J - H)/E(H - K s ) ratio is sen- 
sitive to the adopted shape of the near-infrared extinc- 
tion curve, often tak e n to be a power-law in w avelength 
(|Cardelli et al.l 119891 : iMartin fc Whittetl I1990T ). and to 
the filter bandpasses and the intrinsic spectra of the 
stars being measured, because of color corrections which 
change with increasing extinction. Simulating all of these 
effects, we find that to explain the observed color ex- 
cess ratio, the power-law index for near-infrared extinc- 
tion is 1.85 ± 0.10, encouragingly close to that found to 
be common from studie s of individual re ddened stars 
(|Martin fc Whittetl [1991 . iHe et al.l (fl995T) found a ra- 
tio of 1.64 ± 0.2(23 for highly obscured OB sta r s with 
E{J - K s ) up to about 0.7. Ilndebetouw et al.l (|2005l ) 
obtained a value of 1.73 ± 0.2 using various analyses of 

20 Dobashi also combined his two color excess maps into a rep- 
resentation of Ay ■ For the Vela region at least, we found that this 
is not well correlated with the underlying color excess maps. 

21 Photometry was in the SAAO system but the ratio of 
the color excesses should be similar after color transformations 
to the 2MASS system; see summary by J. M. Carpenter at 
www.astro.caltech.edu / jmc / 2mass / v3 /transformations / 



() 



Martin et al. 



2MASS data probing to the much greater column densi- 
ties accessible in the infrared (E(J—K a ) up to about 2.5). 
The region that we have studied is therefore fairly nor- 
mal, reinforced by the relatively moderate column densi- 
ties in the rectangles chosen (Table [5]) . 

Often near-infrared color excess maps are presented in 
terms of the more familiar measure of extinction Ay , via 
a "total to selective extinction" conversion like 

A v = R XY E(X -Y). (6) 

From studies of individual stars sampling the local 
diffuse interstellar med ium the scaling coefficie nts fre- 
quently adopted, from iRieke fc Lebofskvl (|1985f ) ignor- 
ing the slightly different filter sets, are Rhk = 15.87, 
R JH = 9.35, and R JK = 5.89E3 Even if the shape of the 
near-infrared extinction curve were fairly universal, there 
are significant changes in the shape of the visual to ul- 
traviolet extinction curve, often parameterized in terms 
of Ry = Rbv ; the r atio of total to selective extinction 
(|Cardelli et al.l fl989l) . Increasing Ry from 3.1, charac- 
teristic of the diffuse interstellar medium, to a value of 
5.5, that might be more typical of dark clouds, lowers the 
scaling coefficients by 20%. Furthermore, there are color 
corrections at high column densities. We have no direct 
evidence what scaling to Ay would be valid for these 
Galactic lines of sight toward Vela. Therefore, we decided 
to use the sum of our two maps, E(J — K s ), as the best 
measure of column density, staying close to the directly 
observable color excesses, mitigating any hidden effects 
of grain evolution, and more generally avoiding unneces- 
sary, often hidden, assumptions. Of course, E(J — K s ) 
still has to be calibrated to give Nh (Section [6]) . 

As summarized in Table [1] within the various rect- 
angles measured, E(J — K s ) ranges typically from 0.24 
to 1.1 mag, reaching as low as 0.17 mag in the most 
diffuse field (rectangle 11) and as high as 3.7 mag to- 
ward the main dense cloud in Vela C (rectangle B). In 
a more familiar metric, these color excesses would cor- 
respond to Ay — 0.94, 1.4, 6.5, and 22 mag, using the 
above-mentioned scaling coefficients. However, such scal- 
ing needs to regarded with caution, because the high val- 
ues of color excess and extinction are much beyond those 
toward stars for which the full infrared to ultraviolet ex- 
tinction curve, the underlying dust size distribution, and 
the relationship of Ay to color excesses and TVh have 
been studied directly. 

5. EMISSION AND COLOR EXCESS 

Figure @] illustrates the correlation of /250S the contin- 
uum emission at 250 /im, with the E(J—K S ) color excess. 
Recalling that the dust emission probes the whole line of 
sight, while the depth to which color excess probes is in 
principle limited by both attenuation and sensitivity, the 
correlation is remarkable. The slope of the correlation 
characterizes the emissivity e c (1200) of the dust that is 
causing the spatial variations in column density across 
the map. The zero point of the BLAST maps is not 
needed, but by the same token any fairly uniform screen 
of material cannot be characterized by this correlation 

22 From these scaling coefficients, the ratio of E(J — H) to 
E(H — K B ) is 1.7, and the associated power-law index is 1.61 
HCardelli et al.lfT989T) . 

23 Actually I v but labelled with the wavelength of the passband. 
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Figure 4. Pixel-pixel correlation of 250 /im BLAST emission with 
the color excess E( J — K B ), in rectangle 7. The solid line has a 
best-fit slope e c (1200) = 249.0 ± 5.0 MJy sr _1 mag -1 . The A v 
scale along the top is only approximate, but is provided for context. 

technique. 

The dispersion about the fit in Figure 2] is 0.06 mag 
if measured horizontally or 16 MJy sr" 1 vertically. We 
have estimated the error in the E{ J — K s ) map using 
the scatter about the correlation between the two maps 
E(J — H) and E{H — K s ), examined rectangle by rect- 
angle. From this approach, the rms in the E(J — K s ) 
map is ~ 0.03 mag. Possibly this is an underestimate 
because of correlated errors. Therefore, as an alternative 
we looked at the rms in the adopted E{J—K S ) map in rel- 
atively smooth regions of low color excess and found 0.05 
mag; this could overestimate the error because a part is 
due to real cirrus fluctuations. For comparison with the 
rms for 7 2 50 below, 0.03 - 0.05 mag corresponds to 7 - 
12 MJy sr" 1 . Using similar strategies we have estimated 
the error in /250 : from the scatter about the correlation 
of the J250 map with the /350 map, rectangle by rect- 
angle, we estimate the rms in J250 to be 2 MJy sr _1 ; 
and the rms fluctuation in the dimmest part of the 7 2 50 
map is 3 MJy sr _1 . Returning to Figure HI our estimates 
suggest that the error in E(J — K s ) makes the larger 
contribution of the two and that taken together there is 
no compelling reason to require some additional cosmic 
scatter about the correlation. 

The correlations with the other BLAST bands are 
equally good, but, as foreshadowed by the above cross- 
band comparisons fSection l3~Tj) . the degree of correlation 
of 100 /im emission with E(J — K s ) is less, again possibly 
due to changes in temperature within the volume probed. 

Even though we chose rectangles to avoid strong 
sources, we explored the possibility that clumpy dust 
that would contribute to the BLAST emission map 
might be missed in the sampling of stars underlying the 
E(J — K s ) map. We found no evidence for "excess" emis- 
sion at high color excess even in rectangle B. In the lat- 
ter rectangle there is possibly a slight deficit, related to 



Submillimeter Dust Opacity 



7 



00 
S3 



100 




1000 



100 



Wavelength [//m] 



Figure 5. SED from the slopes of the correlations of dust emission 
with near-infrared color excess E(J — K s ) for rectangle 7. The 
best-fit temperature, 15.2 ± 0.2 K, is similar to that from the 
cross-band correlations (Figure [3j - The SED amplitude e c (1200) 
= 230.3 ± 7.2 MJy sr -1 mag -1 is close to the value obtained from 
direct correlation (Figure 0). 

clumps of cool dust; in that extreme case our assumption 
of a constant temperature for dust within the rectangle, 
rather than systematics in estimating the color excess, is 
probably the limiting approximation. 

5.1. SED 

Fitting the slopes of the correlations between dust 
emission and color excess to a parameterized single- 
temperature SED yields the dust temperature recorded 
in Table [5J An example is given in Figure \5\ These tem- 
peratures from the SED fits are around 15 K, definitely 
cooler than the 17.9 K typical of the high latitude diffus e 
interstellar medium (jPlanck Collaboration et al.ll2011al ). 
In Table [2] we also present the emissivity e c (1200) = 
l25o/E(J—K s ), the best-fit amplitudes obtained from the 
SEDs. Ultimately (Section [7J we will derive the opacity 
from these observables. 

As a consistency check we verified that T derived 
from the l v - E(J — K s ) correlations here is close to 
that obtained via the relative SED from the cross-band 
emission-map correlations (Section |3.1[) . This supports 
our premise that relative changes in column density are 
equally well sampled by changes in BLAST and IRIS 
emission and in near-infrared excess. 

5.1.1. SED weighting and parameter errors 

The slopes of the correlations of BLAST and IRAS 
with E(J — -Ks), and their uncertainties, were obtained 
using SIXLIN. For a given rectangle the uncertainties 
were similar for all BLAST bands and about half of that 
for the IRAS band. This ratio persisted for all rectan- 
gles and so, as in Section 13.11 we have adopted this ra- 
tio uniformly in performing weighted SED fits. We con- 
sider the systematic effect of calibration errors between 
BLAST and IRAS in Section E3 Here a reduced x 2 of 



about unity for the fit typically requires an increase of the 
nominal SIXLIN uncertainties to about 6% and 12% for 
the BLAST and IRAS bands, respectively. This increase 
results in more conservative errors on the model param- 
eters T and e c (1200) (ultimately scaled to <7 e (1200)) and 
on the dependent quantity P (Eq uation ([ 5] )), as calcu - 
lated by Monte Carlo simulation (Cha pin et al.l l2008jh 
A set of 500 realizations of mock data is generated start- 
ing with the actual slopes and adding Gaussian noise 
with the above uncertainties. For each realization the 
SED was fit and the corresponding parameters recorded. 
Finally, the uncertainty on each quantity was obtained 
by fitting a Gaussian to the histogram of the generated 
distribution. These are the statistical errors reported in 
Table [2] By keeping a record of each fit we also tracked 
the correlations of the errors and so can produce the el- 
liptical 1-ct confidence intervals in, for example, the T - 
er o (1200) plane (see Section [7731 . 

6. HYDROGEN COLUMN DENSITY AND COLOR EXCESS 

To calculate the opacity, <j c {v), we need to relate the 
color excess E(J—K S ) to the column density of hydrogen, 
Nr. Often this is obtained by converting E(J — K s ) 
to Ay and then using the ratio Ah /Ay, all for values 
found in the local diffuse interstellar medium. This is 
arguably not justified for the high column density lines 
of sight where it is applied, and therefore is worthy of 
some reflection. 

./Vh toward individual stars has been measured 
using Lyman a absorption and the absorption 
lines of molecular hydrogen using Copernicus and 
the Far Ultraviol et Spectroscop ic E xplorer (FUSE) 
(IBohlin et al.lll978t ISavage et ali l l 977t IDiplas fc Savage! 
119941: iRachford et all 120021 l2009h . But sensitivity re- 
quirements for these ultraviolet measurements have lim- 
ited the column density probed to Ay about 3 mag or 
E(J — K s ) about 0.5 mag, lower than in the fields stud- 
ied here, despite these fields being selected for relatively 
low column density within the Vela map. If the higher 
column density here were simply the result of a long line 
of sight through diffuse material, then it might be ar- 
gued that the material would be similar to what has 
been probed directly using individual stars. However, 
the lines of sigh t within the rectangles are at least par- 
tially molecular (|Yamaguchi et al.lll999f) , and so the ma- 
terial is probably more localized along the line of sight 
with higher spatial density. 

Relating Nh to E(B — V), these studies of individual 
stars have found that 



5.8 x 10 2i E(B - V) cm" 



(7) 



Equation (J7J can be recast in terms of Ay if Ry is known 
for individual lines of sight, but we think that this is 
even less appropriate for our application. Evidence for 
dust evolution on lines of sight passing through dense 
material comes fr om changes in t he op tical-ultraviolet 
extinction curve (Kim & Martin] 119961 and references 
ther ein), which can be p arameterized by the changes in 
Ry ([Cardelli et al.|[l989f) . Because of this, one might ex- 
pect deviations from a simple linear relationship between 
Ah and Ay or even E(B — V), especially for dense re- 
gions with high column density. There is an indication 
of an increase in the slope in the Equation ([7J correla- 
tion to 6.6 x 10 21 cm -2 mag -1 for some higher column 
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Table 2 

Dust emission in the Galactic plane toward Vela 



ID 


T 


e c (1200) a 


rKQ 


<r e (1200) 


P 


E(J — K B ) range 




K 


MJy sr~* mag -1 


2 —1 

cm gm 


10 -25 cm 2 H -1 


10- 31 W H _1 


mag 


1 


16. 8± 0.4 


202. 6± 7.6 


0.09±0.01 


2.0±0.2 


5.3± 0.4 


0.25/0.68 


2 


15. 5± 0.2 


252.1± 6.6 


0.15±0.01 


3.4±0.2 


5.6± 0.2 


0.34/0.85 


3 


15. 2± 0.5 


190. 1± 11.0 


0.12±0.01 


2.8±0.3 


4.0± 0.4 


0.17/0.59 


A 
■I 




OCC Q_l_ C 1 


U. loztU.Ul 


^ n_i_n i 


D.ozt U.z 


n o/i /n 71 

U.Z4/ U. ( 1 


5 


13. 8± 0.7 


83.8± 8.6 


0.08±0.01 


1.8±0.3 


1.5± 0.3 


0.31/0.59 


6 


14.9± 0.2 


164.5± 4.1 


O.lliO.01 


2.6±0.1 


3.4± 0.1 


0.34/0.85 


7 


15. 2± 0.2 


230. 3± 7.2 


0.14±0.01 


3.3±0.2 


4.9± 0.3 


0.51/1.10 


8 


15. 1± 0.2 


181. 6± 5.4 


0.12±0.01 


2.7±0.2 


3.8± 0.2 


0.42/0.85 


9 


15. 1± 0.2 


249.0± 7.2 


0.16±0.01 


3.7±0.2 


5.3± 0.2 


0.25/0.59 


10 


15. 2± 0.4 


124.7± 6.6 


0.08±0.01 


1.8±0.2 


2.6± 0.2 


0.34/0.68 


11 


15. 2± 0.1 


263.6± 5.0 


0.16±0.01 


3.9±0.1 


5.6± 0.2 


0.17/0.42 


12 


15. 1± 0.2 


189.0± 4.8 


0.12±0.01 


2.8±0.2 


4.0± 0.2 


0.34/0.81 


A 


12.4± 0.2 


75.3± 2.9 


0.11±0.01 


2.6±0.2 


1.2± 0.1 


0.42/1.70 


B 


12. 2± 0.3 


99. 9± 5.3 


0.16±0.02 


3.8±0.4 


1.5± 0.1 


1.19/3.74 



a Best fit amplitude from fitting SED 

density lines of sight (Ay of 1 to 5 mag) studied with 
FUSE (|Rachford et al.ll200l [2009h . However, it is not 
known whether this trend is maintained in more dense 
regions where probing the total hydrogen column density 
directly is not possible. 

Because the ratios of near-infrared color excess to JVh 
probably change less significantly as the gra ins evolve 
(|Martin fc Whiketl[l99a Kim fc Martinlll996l) . it seems 
advantageous to examine directly the correlation be- 
tween JVh and our observable E(J — K s ). 

We used atomic as well as molecular h y drogen col- 
umn densities measured b y I Savage et al.l (|1977f ) and 
iRachford et ail (|2002l [20091 ) . B elow a thresh o ld of J V H = 
0.6 x 10 21 cm" 2 , as noted bv iSavage et al.1 (|1977| ). the 
majority of hydrogen is in atomic form. For higher col- 
umn densities, with a mixture of conditions along the line 
of sight, our best- fit approximation to the trend in the 
data above the threshold where both forms of hydrogen 
are measured is 

JV H = 1.52 x (JV H i - 0.6 x 10 21 ) + 0.6 x 10 21 cm" 2 . (8) 

iDiplas fe Savage! (|1994l) measured atomic hydrogen along 
many more lines of sight for which E(J — K s ) can be ob- 
tained, but not molecular hydrogen, and so Equation ((8j) 
was used to make a correction. At any column density 
above the threshold there is dispersion in the fractional 
amount of hydrogen in molecular form and so this cor- 
rection is valid only statisti cally. 

For the program stars of IDiplas fc Savage] |l994), we 
extracted archival 2MASS photometric measurements 
via GA TOFF^. We concentrated on the J and K bands 
to maximize the differential extinction and to take ad- 
vantage of the somewhat better photometry than for the 
H band. There is a bimodal distribution in the reported 
photometric uncertainty, the higher peak relating to sat- 
uration for the brighter O and B stars. Based on the 
lower peaks at J and K, characterizing the normal pho- 
tometric error, we have selected those sources which have 
a combined error in J—K less than ?><jjk , i-e., 0.025 mag. 
To find the color excess E(J — K s ) we used the depen- 
dence of intrinsic colors on spectral classification given 

24 http:/ /irsa. ipac.caltech.edu/applications/Gator/ 
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Figure 6. Correlation between JVh and E( J — K B ) . Solid green 
and r ed circ les represent Nn mea s ureme nts from ISavage et al.l 
119771) and IRachford et~aTl 120021 . I2009T) , respectively. Plus 
symbols are for Nn from atomic hydrogen measurements by 
iDiplas fc Savage (1994), corrected statistically using Equation J5J. 
The best-fit line (solid) has slope (11.5 ±0.5) X 10 21 cm -2 mag -1 . 
The dashed line with slope 11.0 X 10 21 cm -2 mag -1 is from Equa- 
tion (0 converted using standard color ratios (see text). Open 
circles mark the data excluded in the process of iterative fitting. 
The Ay scale along the top is only approximate, but illustrates the 
limited range covered in this calibration compared to that found 
in the color excess map in Figure [21 even within the rectangles 
analyzed. 



bv lStraizvs fc Lazauskaitel (|2009D . 

A preliminary plot oi E(B — V) vs. E(J — K s ) showed 
that some of these selected program stars are much red- 
der in E(J—K S ) than could be expected from interstellar 
extinction. We confirmed from the spectral classifica- 
tion that most of the anomalous stars are known Be or 
emission-line stars; these have near-infrared emission in 
addition to that from the photosphere. Hence, to refine 
our source selection we excluded sources in the color ex- 
cess plane lying beyond 2.5cr from the correlation line 
E(B — V)/E(J — K s ) = 1.9; this precaution serves to 
exclude all of the Be stars from our final list, without 
unnecessarily biasing our results below. 
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Figure [5] shows a good correlation between ATjj and 
E(J — K s ), We have used an iterative fitting technique 
to identify and exclude a few outliers more than 3<r from 
the correlation. Not surprisingly, an d justifying this 
appro ach, the outliers are stars from IDiplas fc Savage! 
(1994), where we have estimated the total hydrogen col- 
umn density from the atomic hydrogen column density 
measurements using Equation ((5J. 

The best-fit line is 



10 21 cm" 2 



(11.5±0.5)£7(J-if s )-0.07±0.01. (9) 



For comparison, conversion of the slope in Equa- 
tion J7J) using standard color ratios appropriate to 
an interstellar extinction curve with Ry = 3.1 
(|Rieke fc Lebofskvi U985t ICardelli et all fl989h gives a 
slope 11.0 x 10 21 cm' 2 mag -1 . This is close to what 
we have found directly, perhaps reassuringly, but also to 
some degree coincidentally, given the range of conditions 
sampled (the converted value would be only half as large 
for Ry — 5.5). Note the significant dispersion about the 
average relation. 

However, it has only been possible to obtain a direct 
calibration of this relation to E{ J — K s ) about 0.4 and 
so applications at higher column densities need to be 
viewed with caution. This includes the present appli- 
cation, where the top of the E(J — K s ) range in most 
rectangles (see Table [2]) is beyond the calibrated range. 
On the one hand, to the extent that the larger values 
are simply the result of longer pathlengths, the calibra- 
tion should stand. But, on the other hand, if the larger 
values of E{ J — K s ) are the result of increased volume 
density then the grains might evolve. For example, cal- 
culations of time-dependent extinction curves resulting 
from gr ain evolution by ic e-mantle formation and aggre- 
gation ()Ormel et al.ll20lTl ) show how E(J — K s ) might 
increase for a given column of material, at least initially, 
which would decrease the ratio Nh/E(J — K s ). This 
would in turn raise the derived opacity (Equation ([ID])), 
exaggerating the changes reported below. 

The same caution about lack of direct calibration also 
holds for any application which uses such measures of 
infrared color excess to gauge column density. The situ- 
ation is further muddied, unnecessarily, when the column 
density is cited in terms of Ay. 

7. RESULTS 

7.1. Opacity at 1200 GHz or 250 [im 

Using the temperatures and the amplitudes from the 
SED fits in Section[S] together with the ratio N^/E(J — 
K s ) from Section [51 we calculated the opacity from 

a c (1200) = e c (1200)/[B„ o (T) x N K /E( J - K s )}. (10) 

Recall that er o (1200) = (rninrno. The derived values 
are recorded in Table [5] along with their Monte Carlo 
errors (Section I5.1.1[) . The typical opacity in these re- 
gions is about 2.8 x 10 -25 cm 2 H _1 or equivalently 
0.12 cm 2 gm" 1 . There are considerable variations above 
what can be accounted for by the errors. Furthermore, all 
values are significantly above what is typical of the high 
latitude diffuse interstellar medium , 1.0 x 10 -25 cm 2 H _1 
(Planck Collabor ation et aL1l2011al ). 



Rectangle A coinc ides with cloud 28 of 
lYamaguchi et al.l (|1999l ) for which iNetterfield et al.l 
(|2009D have estimated the dust opacity to be 
rno — 0.16 cm 2 gm" 1 by comparing the integrated 
submillimeter dust emission with the total mass of gas 
estimated from the CO emission. The latter introduces 
an uncertainty of a factor 1.5 - 2. Our new value is 
0.11 ±0.01 cm 2 gm" 1 . 

7.2. Integrated emission 

A physical quantity of interest is the total energy emit- 
ted by dust per hydrogen atom, P (Equation (|SJ)). For 
the diffuse high Galactic latitude interstellar medium 
the value is fairly uniform near 3.8 x 10~ 31 W H _1 
(jPlanck Collaboration et aL1l2011al) . The typical value 
found here (see Table [3]) is somewhat higher, 4.5 x 
10 -31 W H _1 . However, not surprisingly, there is consid- 
erable variation in the Galactic plane (see also Section lST2l 
and Figure [5] below) . 

7.3. Relationships 

The parameters T, ct o (1200), and P for any rectangle 
are related at a fundamental level through Equation ([5J . 
Because we have used a fixed /? = 1.8, this relationship 
can be quantified as 



(T/T ) 5 - 8 



(P/P ) (a e (1200)/a ) 



(11) 



using the above-mentioned high latitude diffuse ISM val- 
ues for normalization (Section [5]). 

The parameters derived for the different rectangles 
(Table [2]) and their elliptical 1-er confidence intervals 
(Section I5.1.1|) are displayed in two complementary di- 
agrams, T - cr e (1200) (Figured]) and P - a c (1200) (Fig- 
ure |5]). These figures include loci according to Equa- 
tion (1111) along which the third parameter is constant. 
Because of fixed /3 and Equation (|TT]) , the third possible 
diagram, T - P, contains no independent new informa- 
tion. 

7.4. Other errors 

Throughout this study we have avoided high extinction 
regions, owing to a larger uncertainty in the color excess 
and in its conversion to A^h- 

The adopted ratio of Nb_/E(J — K s ) might not be uni- 
versally applicable. A systematic 5% uncertainty in this 
ratio, or a larger error for individual measurements given 
the dispersion about the relation (Figure [6] also relevant 
to Equation ([7]) which is just the average), induces (in- 
versely) errors of the same size in cr e (1200). This is poten- 
tially a major contribution to the error budget. While 
we might have minimized this uncertainty successfully 
by avoiding high extinction regions, by the same token 
our derived values of opacity might not be representa- 
tive of the unanalyzed high extinction regions, including 
compact sources. Because of the effects of grain evolu- 
tion mentioned in Section [6] (see also Section I&3]) . any 
application where near-infrared color excess is used to 
estimate high column densities needs to be assessed crit- 
ically. We have investigated whether part of the trends 
seen in Figures [7J and [5] might be induced by changes in 
fj c (1200) originating in a trend of N-r/E(J — K s ) with 
environment; however, a scatter plot of <r c (1200) vs. the 
average < E(J~K S ) > for the different rectangles shows 
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Figure 7. Dust temperature, T, vs. emission opacity, <r c (1200). 
Values for the regions studied here in the Galactic plane toward 
Vela are represented by the red filled squares with elliptical 1-cr 
confidence inte rvals. Loci are for constant P, the third parameter 
(Equation <TTt): solid 3.8 X 10~ 31 W H" 1 for the diffuse high lat- 
itude ISM (Planck Collaboration ct al. 20fla); and dashed, three 
times less. If /3 were fixed at 2.0 in the SED fit, each solution shifts 
systematically at roughly constant P; this is illustrated for rect- 
angle 7 by a line connecting to the star. Also for rectangle 7, the 
systematic effect of a 15% calibration error between BLAST and 
IRIS is illustrated by the dotted line segments (see text). The filled 
circle locates the standard values fo r the diffuse high latitude ISM 
(Planck Collaboration ct al. 2011a). Values for individual diffuse 
high latitude regions (Planck Collaboration ct al. 2011a) are also 
plotted, with black and blue distinguishing between LVC and IVC 
components, respectively. Some further comparisons with empir- 
ical results (Appendix [B} are plotted. Open ci rcles show Planck 
results from Planck Collaboration ct al. (2011b) toward the Tau- 
rus molecular cloud for lines of sight both atomic (upper left) and 
molecular; the vertical error bars indicate the ranges of tempera- 
ture. The fille d diamond is anothe r estimate for the Taurus region 
using Spitzer (Tcrcbcv et al. 2009). The squares show the range of 
values found using Hersch el data in the environs of a Planck cold 
clump ( Juvela ct al. 2011); these results support the trend with T 
(and column density) established by the other data and extend it 
to even higher <j e (1200) at the highest column densities. 
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Figure 8. Like Figure[7]but for power emitted by dust per hydro- 
gen a tom , P, vs. emission opacity, <r c (1200), and with loci (Equa- 
tion (TTJ) for different T: 13, 15, 16, and 17.9 K from lower to 
upper. 



no definitive trend, as can be judged by inspection of Ta- 
ble [2] as well. 

We kept the rectangles fairly small so that the assump- 
tion of a uniform temperature for the dust was at least 
plausible; however, this does not control for temperature 
changes along the line of sight. 

There are calibration errors which could produce sys- 
tematic effects. For the BLAST bands these are 8.1, 
7.1, and 7.8% at 250, 350, and 500 fj,m, respectively, 
and by the nature of the calibration technique they ar e 
well correlated between the bands (|Truch et al.l [2009) . 
The c alibration error of the IRIS 100 Mm image used is 
13% ()Miville-Deschenes & Lagachdl2005l ). Because the 
calibration techniques are independent, there is the pos- 
sibility of a systematic effect when SEDs are fit to the 
combined data. If IRIS were adjusted upward, then T 
would be higher and <r e (1200) lower. If instead it were 
BLAST that was adjusted downward, by the same rel- 
ative amount, then the same higher T would be found 
and <T (1200) would be even lower. Such a two-part tra- 
jectory is illustrated in the figures for rectangle 7, for a 
relative calibration error lowering BLAST/IRIS by 15%. 

For consistency and comparison to other studies we 
adopted (3 = 1.8. The choice of j3 affects the parame- 
ters derived from the SED fit in a systematic way. The 
total power radiated P is quite insensitive to changes in 
/3, because the SED is required to fit the data across 
the whole range being integrated. However, both the 
amplitude and temperature change, the latter affecting 
the derived opacity more profoundly because of the non- 
linearity of the Planck function (Equation (fTU|)'). For 
example, if (3 were increased from 1.8 to 2.0, the tem- 
perature would be systematically lower and the opacity 
systematically higher. This shift is illustrated in Figure[7] 
for rectangle 7. Such effects seem unlikely to account for 
the much larger values of the opacity found here com- 
pared to those in the high latitude interstellar medium. 
On the contrary, in Appendix lAl we explore the potential 
effect of /3 further, in particular through a possible /3 - 
T relationship, and conclude that this would exaggerate 
the changes in opacity already found with fixed /?. From 
this perspective, the estimated magnitude of the changes 
found and discussed below is conservative. 

8. DISCUSSION 

There are significant variations from rectangle to rect- 
angle in the derived parameters T, er c (1200), and P. 

For comparison with these new results for the Galactic 
plane toward Vela, the filled circle in Figures [7] and [5] lo- 
cates the stand ard values for the diffuse high Galact ic lat- 
itude ISM from lPlanck Collaboration et all (|2011aD that 
were used to normalize Equation (fTTj) . The data plotted 
for individual diffuse high latitude regions are from SED 
fits (0 = 1.8) to the emissiyities r eported in Table 2 of 
iPlanck Co llaborati on et ail ()2011aD . in which local veloc- 
ity clouds (LVC) and intermediate velocity clouds (IVC) 
have been separated. Note that the dust emission un- 
derlying these SEDs is measured from a combination of 
Planck and IRIS data. Even at high latitudes there is 
considerable variation, but interestingly the ranges of T 
and er c (1200) do not overlap with those found in the 
present study. 

Overall there is a trend in Figure [7] of decreasing T 
with increasing <t c (1200). Although the error ellipses are 
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aligned roughly along this trend, the values span a much 
larger range than can be attributed to the individual er- 
rors of the independent measurements. The range of P 
in the present study extends to both higher and lower 
values than seen at high latitudes. 

Other estimates discussed in Appendix |5] are plotted 
in the figures. These support the trend of decreasing T 
with increasing er c (1200) seen above. 

It is useful to think of Equation (jTTJ) in terms of cause 
(on the right hand side) and effect, T. The P calcu- 
lated here is for emission by big grains in thermal equi- 
librium and is therefore also the total power absorbed by 
these grains when exposed to the interstellar radiation 
field (ISRF). The opacity cr o (1200) measures an intrinsic 
property of the big grains, how efficiently they can emit. 
This emission opacity and the absorbed P that needs 
to be emitted determine what the equilibrium T must 
be. Thus it can be seen that big grains will be cooler 
in a less intense ISRF, and/or if they were to evolve to 
have a lower absorption opacity and/or a larger emission 
opacity. 

Here we examine the evidence for such changes using 
the complementary diagnostic Figures [7] and [8] Consid- 
ering T as an effect, we concentrate further discussion on 
the causes, er c (1200) and P. 

8.1. Changes in er o (1200) 

A main goal of this paper was to quantify the dust 
emission opacity in a new environment near the Galac- 
tic plane, that is of higher column density than the high 
latitude ISM, and at least partly molecular. We have 
found that in this environment cr c (1200) is higher by 
typically a factor three and, exte nding the finding by 
iPlanck Collaboration et all (|2011af ). it changes from re- 
gion to region. 

Changes in the emission opacity are certainly intrigu- 
ing but not understood. One possibility is that the opac- 
ity is raised when grains aggregate, changing the basic 
structure to something more porous and fractal than ho- 
mogeneous ( Oss enkopf fc Henningl [1994 ). This has been 
discussed for dense molecular clouds, where grains also 
develop ice mantles ()Ormel et al.ll201l[) , but its relevance 
to such evolutionary changes occurring even within the 
more diffuse medium seems less obvious. Perhaps one 
needs a change in perspective on the direction of grain 
evolution, regarding the higher values found here as "nor- 
mal" for dense regions and the lower values as the result 
of evolution of dust back toward a diffe r ent st ate in the 
diffuse interstellar medium. See Uonesl (|2009D for a re- 
lated discussion on extinction curves. 

8.2. Changes in P 

One of the possible reasons for the range of values of 
P is variation of the interstellar radiation field (ISRF) 
in the Galactic plane. Attenuation in dense molecular 
clouds seems the likely cause of the lower P in rectangles 
A and B. Likewise, it is at least plausible that the ISRF 
is higher in those regions (rectangles 2, 4, 9, 11) with 
P significantly larger than 5 x 10 -31 W H _1 . However, 
as can be seen from Table [2j a scatter plot of P vs. the 
average < E(J — K s ) > for the different rectangles does 
not show any definitive trend. 

Especially for rectangles A and B, with the most ex- 
treme conditions, we have to be aware of other additional 



uncertainties, such as N-r/ E(J — K s ) being different than 
adopted; a decrease in this ratio (Section [6J would raise 
both the derived a (1200) and P (deduced from the emis- 
sion) proportionately, at constant T (lines in Figure [8]). 

Another factor is grain evolution. If grains are evolving 
enough to change the emission opacity cr o (1200) signifi- 
cantly, it is at least plausible that the capacity to absorb 
is also changing. However, there are no near-infrared to 
ultraviolet extinction curves from which to quantify such 
changes. The evolutionary effects are important to un- 
derstand because the absorption opacity directly affects 
the power absorbed from the ISRF and hence P observed 
in emission. For the typical sizes of big grains in the 
diffuse interstellar medium, the absorption cross section 
is approaching the geometric cross section. With grain 
growth by accretion and aggregation, this ratio of cross 
sections saturates at a value near unity. Grain growth 
also increases the mass faster than the geometric cross 
section, driving down the absorption opacity. Modeling 
these competing effects would take a detailed grain model 
and theory of grain evolution, as well as an accounting of 
the spectral shape of the (attenuated) ISRF, well beyond 
the scope of this paper. 

Another potential factor is the dust-to-gas mass ra- 
tio. All of the fields with data in these figures are at 
the solar Galactocentric distance, and so the underlying 
metallicity is likely the same. Furthermore, the depletion 
is already high in the diffuse ISM, leaving little room for 
dust mass to increase in more dense regions; however, ice 
signatures do appear in molecular clouds. On the other 
hand, there is independent evidence for reduced deple- 
tion in some IVCs, which would lower both cr e (1200) and 
P at constant T. There is a hint of such an effect in 
Figured 

8.3. Insight from theoretical modeling of grain evolution 

The emission cross section of dust grains in high lati- 
tude diffuse interstellar clouds is better constrained em- 
pirically than in translucent (Ay in the range of 1 to 5 
mag) or dense molecular clouds ( Appendix IE]) . 

In the ir comprehensive revi ew of existing empirical es- 
timates, [Henning^tjil] (|1995f) comment that "This large 
scatter [in submillimeter opacities] is very probably not 
only related to problems with the observational deter- 
mination of the opacities but may also reflect real dif- 
ferences of the dust populations in different environ- 
ments caused by evolutionary effects (e.g., coagulation 
or accretion of mantles)." This has motivated model- 
ing of the grain evolution and the attendant changes in 
opacity. After reviewing their theoreti cal calculations 
relati ng to the evolution (ice m antles: iPreibisch et al. 
1 19931 : mantles and coagulat ion: lOssenko pf & Hcnning 
1994), iHenning et al.l (|1995| ) recommended values for 
three notional environments with the opacity increas- 
ing in magnitude fr om the diffuse interstellar medium 
(|Draine &; Ledfl984t like current estimat es) to protostel- 
lar c loud envelopes (nn ~ 10 5 cm -3 , IPreibisch et al.l 
Il993t like the values found here for less dense molec- 
ular regions) and even further in dense and cold cores 
(nn > 10 7 cm" 3 , lOssenkopf fc Henninel 119941 : condi- 
tions well beyond what is probed here). Although it 
is acknowledged that the models are instructive rather 
than definitive, the theoretical estimates have often been 
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adopted for the analysis of submillimeter data because 
of the lack of empirically-calibrated opacities for these 
environments; the situation is improving somewhat (Ap- 
pendix[B|) but the most dense and evolved regions remain 
challenging. 

With deeper targeted surveys of individual clouds us- 
ing infrared cameras, some variations of the ratio E(J — 
H)/E{H — K s ) have been found, suggestive of grain evo- 
lution. For examplc F^I the ratio foun d is 1.76 ± 0.07 in 
the p Oph clou d (IKenvon et al.lll998fl, 1.96 ± 0.04 in the 
Cham I cloud (IGomez k. Kenvonll200ll). a nd 2.25 ± 0.07 
in Coalsack Globule 2 (jRacca et al.ll2002f) . In these in- 
vestigations, E(J — K s ) ranged up to 5.65, 3.04, and 
2.44, respectively, all well beyond the top values in the 
rectangles considered here. As emphasized in Section 
Nn/E(J — K s ) has not been calibrated directly for such 
high column densities. 

The multi- waveleng th complexity is hig hlighted by the 
recent calculations bv lOrmel et al.l (|2011l ). motivated by 
evidence for changes in the near-infrared extinction for 
inferred column densities up to iVn ~ 6 x 10 22 cm -2 (no- 
tionally Ay ~ 32 mag), again well beyond that probed 
here. Following the effects of ice-mantle formation and 
(subsequent) grain coagulation, they model the opac- 
ity across the whole spectrum from the ultraviolet to 
submillimeter. As with the earlier models, these results 
show how the evolution can produce not only increases 
in submillimeter opacity but also accompanying dramatic 
changes in the near-infrared and ultraviolet opacity. 

The combined effects of magnitude and slope changes 
in the near-infrared opacity can alter N-r/ E(J —K s ). De- 
pending on the model, significant near-infrared changes 
might develop even before a substantial change in sub- 
millimeter opacity. In the Vela region that we have ana- 
lyzed above there is no change in the infrared slope, and 
assuming no change in the absolute amount of the extinc- 
tion either (Section |6]) we find that there is an increase 
in submillimeter opacity relative to that in the diffuse 
interstellar medium. If the magnitude of the infrared 
opacity, which we cannot measure directly (but see Fig- 
ure , has actually increased, then the derived opacity 
would be even high e r (see Equation (jTUll ). 

The lOrmel et al.l (|2011h results also show a decrease 
in optical-ultraviolet opacity as the grains evolve, a re- 
minder that Ay might not be a good surrogate of column 
density in these evolved regions. This optical-ultraviolet 
opacity decrease would decrease the energy absorbed and 
thus the observed P in emission. However, depending 
on the details of the evolution, its time development, 
and the spectral shape of the radiation field, this de- 
crease might be compensated by an increase in the near- 
infrared opacity. As discussed above, a mis-calibration 
of Nn/E(J - K s ) would scale P and cr o (1200) equally, 
moving the derived quantities along lines of constant T 
in Figure [8] 

9. CONCLUSION 

We have correlated the diffuse interstellar dust emis- 
sion in the Galactic plane toward Vela (BLAST images 
at 250, 350, and 500 /jm and the IRAS image at 100 /im) 
with a map of near-infrared color excess made using 

25 For consistency, we have transformed the original photometry 



2MASS data. Fourteen regions of moderate column den- 
sity were analyzed. The conversion of color excess to col- 
umn density has been examined critically. From stellar 
data we have measured Nn/E(J — K s ) to be 11.5 x 10 21 
cm 2 mag -1 with a considerable dispersion (Figure [6]). 
From the spectral energy distribution of the dust emis- 
sion, we have quantified important properties of the big 
grains, namely the equilibrium temperature T of the big 
grains that are in thermal equilibrium with the inter- 
stellar radiation field (ISRF), their submillimeter opacity 
er c (1200) (the emission cross section per H nucleon), and 
P, the total power radiated per H nucleon. We find that: 

1. The submillimeter opacity is consistently larger than 
for dust in the local high Galactic latitude interstel- 
lar medium, by a factor 2 to 4 relative to the stan- 
dard <7 G (1200) = 1.0 x 10~ 25 cm 2 IT 1 value (re = 
0.043 cm 2 gm _1 ). This is strong evidence for grain 
evolution. 

2. The range of P extends to both higher and lower 
values compared to that found at high latitudes, 
3.8 x 10~ 31 W Br 1 (1.2 L©/M ). This range in part 
reflects variations in the interstellar radiation field. It 
is also influenced by evolutionary changes in the dust 
absorption opacity. In turn, all of the above changes 
lead to changes in the observable equilibrium T . 

3. Compared to the local high latitude dust temperature 
(17.9 K), in this direction in the Galactic plane the 
dust temperatures are significantly colder, typically 
15 K. Somewhat lower temperatures still are found in 
more dense higher column density regions where the 
ISRF is more strongly attenuated. 

4. Continuing the trend found in high latitude fields, 
there is an inverse correlation of T with cx c (1200). 

The recognition that there are changes in the emission 
opacity raises a particular point of caution, because the 
value adopted impacts directly all column densities de- 
duced from dust emission maps, and the masses of com- 
pact structures (clumps, cores, filaments, ridges) within 
them. While values typically being adopted (see Ap- 
pendix |5]) are within the range that we find, we will 
need to understand the underlying causes of the varia- 
tions already observed in order to assess whether there 
are further changes to be expected in the important even 
denser environments where the opacity has not been cal- 
ibrated. 
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Columbia Scientific Balloon Facility (CSBF) staff for 
their outstanding work. Finally, we appreciate the care- 
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which has led to some clarification and elaboration of the 
analysis. 
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APPENDIX 

A. EXPLORATION OF THE IMPACT OF A - T RELATIONSHIP 

There is a considerable lite rature on a possible inverse relationship between (3 and T. However, as is clear from the 
many examples in Figure 3 of Paradi s~et al.l (|2010f) . there is no consensus on the details of this dependence. In practice, 
we do not have sufficient multi- wavelength data to treat f3 as an additional free parameter in the SED fit; without the 
fit being well over-constrained, undue sensitivity could develop to the weighting of the data and issues of calibration, 
for example. Nevertheless, we have gone through the exercise of treating /3 as a free parameter. For our rectangles we 
found that the values of free-/? tended to be a bit larger than 1.8, opposite to the lowering of the apparent f3 expected 
if the SED is broadened because of dust of different temperatures superimposed along t he line of sight. 

Over all of the regions examined here plus those in lPlanck Collab oration et ahl (|201 laf ) there is a considerable range 
of T and a suggestion that f3 is inversely related. This trend could be parameterized as 

/3/1.8 = (T/17.9K)" 2/3 , (Al) 

the power- law index being intermediate among the examples summarized by iParadis et al.l (|2010| ) . 

For the purposes of illustration, following this possible trend through to its potential consequences, we have performed 
SED fits subject to the added constraint of Equation (jAll) . The derived parameters are displayed in Figure O Loci for 
constant P under the same constraint are plotted for reference. The effect of constraining the SED fit with this (3 - T 
relationship is systematic. As expected, each solution remains near the same value of P found for fixed /?. However, 
compared to Figure [71 the values of the other two parameters are spread out over a larger range on either side of the 
fiducial values corresponding to (3 — 1.8 and T = 17.9 K. 

This spreading exaggerates the differences in a c (1200) found from region to region using fixed (3. We are not 
actually persuaded that there is a f3 — T relationship, but in any case we conclude that the estimated magnitude of 
the environmental variations as found for fixed f3 = 1.8 and discussed in the paper is conservative. 

B. VALUES OF OPACITY IN THE LITERATURE 

As in Section [21 we standardize on 1200 GHz (250 jum) as the fiducial reference frequency Vq, a value quite appro- 
priate to on-going analyses of Herschel data, for example. The opacity cr c (1200) and the product rua can be used 
interchangeably through Equation ([J}, as in the lower and upper axes in Figures [7J and HJ the scale factor between the 
respective numbers on the axes is 0.043. Scaling to and from the value of opacity at another frequency depends on (3 
in which there is some uncertainty; we use j3 = 1.8. For example, scaling to 1000 GHz lowers the opacity by a factor 
1.4. 

B.l. Adopted 

A value representative of the iPreibisch et al.l ()1993f ) theoretical core-mantle grain evolution results for conditions in 
the cloud envelopes of prestellar cores, corresponding to 0.1 cm 2 gm _1 at 1000 GHz or rug — 0.14 cm 2 gm" 1 , has 
been adopted by the Herschel G uaranteed Time Key Programs on star formation (Gould Belt Survey, lAndre et al.1 
120101 HOBYS. iMotte et al.l 120 lOh . allowing consistent comparisons between different regions analyzed. Whether this 
is valid for the dense cores extracted is uncertain. 

Analyses based on data (on l y) at lower frequencies, for example by iKerton et al.l (|2001l ) for data from SCUBA 
at JCMT and by IMotte et ail i|2007f ) for data from MAMBO at IRA M, often adopt an opacity a t a lower fiducial 
frequency. These are based on theoretical estimates for hh ~ 10 5 cm~ 3 (Oss enkopf fc Hennindll994l) . and would scale 
compatibly using to /? = 2 to vkq = 0.23 cm 2 gm _1 (or 0.18 cm 2 gm _1 for (3 — 1.8). This is very similar to the value 
adopted for Herschel because of the closely related theoretical basis. 

Analyzi ng the compact sour c es in the Sc i ence Demonstration Phase (SDP) fields of the Open Time Key Program 
Hi-GAL (|Molinari et al.l I2010D . lElia et al.l (|2010[) adopted not only a low fiducial frequency (230 GHz) but also a 
variable j3 in fitting the SEDs; thus comparisons of derived masses are not so straightforward, although effects of f3 on 
T and the implied rno somewhat cancel in the mass estimates. 

In most BLAST papers on Galactic star forming regions (|Chapin et all [200l iTruch et al.l |2008; 
iRivera-Ingraham et al.|[2010t IRov et al.ll2011al [bT) we adopted r« = 0.10 cm 2 gm" 1 ()Hildebrandl 119831: derived em- 
pirically from limited observations of a single molecular cloud, for which the value is "probably good within a factor 
three or four," this has nevertheless often been adopted as a "canonical" value). In the BLAST analysis of the Vela 
molecular ridge region (INetterfield et al.l [200l fOlmi et al.ll2009h we estimated 0.16 cm 2 gm 1 using a CO calibra- 
tion (Section IT. 1[) . For the diffuse interstellar medium toward the CasA supernova remnant (jSibthorpe et al.l 120091 ) 
we used 0.05 cm 2 gm" 1 . A low value like this is also implicit in any analysis with the standard grains in DustEM 
(jCompiegne et alj|201lh . 

B.2. Empirical 

The value of the opacity for the diffuse high latitude interstellar medium is the best determined, through correlations 
of the dust emiss i on wi th Ah from observations of the 21 cm emission line. Using observations from DIRBE and FIRAS, 
iBoulanger et al.l (| 1996T) obtained an emission cross-sectio n a e ( 1200) = 1 x 10~ 25 cm 2 H _1 , and an equili brium temper- 
ature of 17.5 K, in agreement with the value obtained bv lDraine fe Led (|1984l ). Also using COBE data. lLagache et ahl 
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Figure 9. Dust temperature, T, vs. emission opacity, ct o (1200), from SED fits under the constraint of the f3 - T relationship in 

Equation {Al}. Loci are for constant P under the same constraint, for the values used in Figure [7] Compared to Figure the solutions 
for the parameters are spread out over a larger range about the fiducial values corresponding to j3 = 1.8 and To = 17.9 K. However, the 
increased spread is systematic, each point remaining near the same value of P found for fixed ft. Because of this spreading, the coverage of 
this figure has been expanded relative to that of Figure [7] (dotted outline) and even so the most extreme points, an IVC and rectangle B, 
lie off the plot to the upper left and lower right, respectively. For rectangle 7, the star represents the so lutio n obtained when /3 is treated as 
a free parameter in the SED fit; /3 is found to be 2.2, slightly above the trend described by Equation 1AH . and so that solution is similar 
to the constrained result here. 



(1999) found 0.87 ±0.09 x 10~ 25 cm 2 H" 1 , compatible also with the result obtain ed by IWeingartner k, Drainel (|2001f) . 
In the text here we cited and plotted tkq = 0.043 cm 2 gm _1 for T = 17.9 K (|Planck Colla boratio n et "all 1201 lab . 
Recently, multi-wavel ength analysis of correlations of P lanck and IRAS dust maps covering submillimeter to far- 
infrare d wavelengths ([Planck C ollaborati on et al.l 1201 laf ) with new higher-resolution observations of Nh from the 
GBT (IBoothrovd et all 120111 Martin et al. in preparation) showed that there were regional variations about this 
value (see Figure [7]). The atomic hydrogen column density associated with these high latitude clouds ranges over 
Ng ~ 0.2 to 6 x 10 20 cm" 2 which is equivalent to Ay ~ 0.01 to 0.3 mag. For the atomic phase in the Taurus field, 
iPlanck Collaboration et all ()2011bl ) obtained a similar opacity 1.14 ± 0.2 x 10 -25 cm 2 H _1 using Planck data, for 
column densities up to 3 x 10 20 cm" 2 (Ay about 1. 6 mag). 

In their Section 5.4.1, IPlanck Collaboration et all (|2011bD revie w earlier estimates of opacity in more molecular 
regions (see also Section 5.2 in ljuvela et al.l 120111 and Figure 3b in Kramer et al.l 120031 ) . For the molecular phase in 
the same Taurus field they obtained a higher opacity, 2.3 ± 0.3 x 10~ 25 cm 2 H _1 , using Planck data and gauging the 
column density using 21 cm and CO emission-line observations (Nn ~ 1 x 10 22 cm -2 ). Values for both the atomic 
phase and molecular phase are plotted in Figure they follow basically the same trend established by the other data, 
on a locus of slightly lower p ower (see also Figure |SJ . 

Also in the Taurus region, iTerebev et al.l (|2009f) find a similarly higher opacity, by correlating dust optical depth 
(T = 14.2 K) from MIPS imaging from Spitzer, extending to 160 /xm, with the column density. To estimate the 
column density they used near-infrared color-excess (which has the same uncertainties as elaborated here in Section [6]), 
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making a careful analysis of Ay. Transforming their 160-^m opacity to our fiducial frequency, we find cr c (1200) = 
2.1 x 1CT 25 cm 2 H" 1 for (3 = 2 (2.3 for /3 = 1.8); see Figures [JJ and HI iFlagev et all (|2009fl also analyzed similar data 
over a slightly bigger map in Taurus, finding a very similar temperature and opacity. 

With the advent of submillimeter mapping by Planck and Herschel (with zero-point offsets from Planck), SEDs 
have been fit pixel by pixel, resulti ng in maps of dust o ptical depth r and T. From the slope of r versus Nn estimated 
from observations of H I and CO, iBernard et al.l (|2010f ) find values tkq = 0.094 and 0.14 cm 2 gm _1 for the Hi-GAL 
SDP fields at i = 59° and I = 30°, respectively. 

Dividing the r map directly by a map of column density produces a map of opacity (Equation (j4])). Note that for 
this application column density is often obtained by converting a near-infrared color excess (unnecessarily expressed 
as Ay) into Nn and so the very same caution (Sections lr?l and 153]) as to the lack of direct calibration at high column 
density applies, even more so. This pixel by pixel approach has the advantage of tracking opacity changes at higher 
spatial resolution compared to the correlation analyses used here. These spatial changes can be related to changes in 
T as well, e.g., producing a scatter-plot version of Figure [Jj However, it is worth recalling that this method essentially 
assumes that the properties of the dust (T, opacity) are uniform along the line of sight, which might not be the case 
when the properties are ap parently changing si gnificantly in the transverse direction. Mapping the environs of Planck 
cold clumps with Herschel, Uuvela et al.l ([201 1[ ) find opacities typically 0.1 cm 2 gm _1 where T ~ 17 K (f3 = 2), but 
that on the high column density lines of sight T drops to 14.5 K and the opacity rises to 0.2 to 0.3 cm 2 gm . This 
again supports the trend found in Figure [Jj now on a locus of slightly larger power. 
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